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ABSTRACT 



We investigate various problems related to the theory of the electroweak 
phase transition. This includes determination of the nature of the phase tran- 
sition, discussion of the possible role of the higher order radiative corrections 
and the theory of the formation and evolution of bubbles of the new phase. 
We show, in particular, that no dangerous linear terms in the scalar field 
appear in the expression for the effective potential. We have found, that for 
the Higgs boson mass smaller than the masses of W and Z bosons, the phase 
transition is of the first order. However, its strength is approximately 2/3 
times less than what follows from the one-loop approximation. 

The phase transition occurs due to production and expansion of critical 
bubbles. Subcritical bubbles may be important only if the phase transition 
is very weakly first order. A general analytic expression for the probability of 
the bubble formation is obtained, which may be used for study of tunneling 
in a wide class of theories. 

The bubble wall velocity depends on many factors, including the ratio of 
the mean free path of the particles to the thickness of the wall. Thin walls 
in the electroweak theory have a nonrelativistic velocity, whereas thick walls 
may be relativistic. 

Decrease of the cubic term by the factor 2/3 rules out baryogenesis in the 
minimal version of the electroweak theory. Even though we concentrate in 
this paper on the phase transition in this theory, most of our results can be 
applied to more general models as well, where baryogenesis is possible. 



1 Introduction 



Twenty years ago David Kirzhnits discovered that the symmetry between 
weak and electromagnetic interactions should be restored at the very early 
stages of the evolution of the universe Symmetry breaking between weak 
and electromagnetic interactions occurs when the universe cools down to a 
critical temperature Tc ~ 10^ GeV. His results were confirmed by an investi- 
gation performed in 1974 by Weinberg, Dolan and Jackiw and by Kirzhnits 
and Linde , and soon the theory of the electroweak phase transition became 
one of the well established ingredients of modern cosmology. Surprisingly 
enough, a complete theory of this phase transition is still lacking. 

In the first papers on this problem it was assumed that the phase transi- 
tion is of the second order 0, Q]. Later Kirzhnits and Linde showed that 
in the gauge theories with many particles, and especially with particles which 
are much more heavy than the Higgs boson 0, one should take into account 
corrections to the high temperature approximation used in |Q, |^ . These cor- 
rections lead to the occurrence of cubic terms ~ g^(f)^T in the expression for 
the effective potential V{(j),T). As a result, at some temperature, V acquires 
an extra minimum, and the phase transition is first order 0. Such phase 
transitions occur through the formation and subsequent expansion of bubbles 
of the scalar field inside the symmetric phase (f) = 0. A further investiga- 
tion of this question has shown that the phase transitions in grand unified 
theories are always strongly first order [Q. This realization, as well as the 
mechanism of reheating of the universe during the decay of the supercooled 
vacuum state suggested in 0, ^, played an important role in the development 
of the first versions of the inflationary universe scenario P . (For a review of 
the theory of phase transitions and inflationary cosmology see Ref. 0.) 

For a long time it did not seem likely that the electroweak phase transition 
could have any dramatic consequences, unless the Higgs boson is exception- 
ally light. Even though the possibility of a strong baryon number violation 
during the electroweak phase transition was pointed out fifteen years ago by 
Linde [|] and by Dimopoulos and Susskind , only after the groundbreaking 
paper by Kuzmin, Rubakov and Shaposhnikov |Ty] was it realized that such 
processes do actually occur and may erase all previously generated baryon 
asymmetry of the universe. 
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Recently, the possibility that electroweak interactions may not only erase 
but also produce the cosmic baryon asymmetry has led to renewed interest 
in the electroweak phase transition. A number of scenarios have been pro- 
posed for generating the asymmetry [|r^ - |[TB|. All of them require that 



the phase transition should be strongly first order since otherwise the baryon 
asymmetry generated during the phase transition subsequently disappears. 
In all of these scenarios the asymmetry is produced near the walls of the 
bubbles of the scalar field 0. Thus, an understanding of the nature of the 
phase transition in the electroweak theory and investigation of the structure 
of the bubbles produced during the phase transition are of some importance. 
For this purpose one should make a much more thorough analysis of the 
electroweak phase transition than the analysis which is necessary for an ap- 
proximate calculation of the critical temperature. 

Unfortunately, despite the fact that one is dealing with a weakly coupled 
theory, many aspects of the phase transition are surprisingly complicated. 
Indeed, the literature contains contradictory claims and statements on almost 
every important question. In this paper, we will attempt to resolve a number 
of these questions, or at least to delineate the issues which are crucial to a 
complete analysis. We will confine our attention to weakly coupled theories, 
not because strongly coupled theories (e.g. technicolour theories) are not 
of interest, but simply because we will find these problems quite difficult 
even in theories with explicit Higgs particles. Even more specifically, we will 
discuss the simplest version of the electroweak theory containing only one 
Higgs boson. We doubt that the baryon asymmetry can be generated in this 
simple model; in fact, we will find some new evidence against it. However, 
this model will help us to illustrate various possibilities which may be realized 
in more complicated theories. 

The first problem to be studied is whether the phase transition is first 
order, and precisely how strongly first order it may be. One clear requirement 



arises if no net B — L is produced at the transition (in the model of |T6[ a net 
-B — L is produced). In this circumstance, as first stressed in []TT|, the rate of 
baryon number violating transitions after the phase transition is completed 
must be smaller than the expansion rate. In practice, this means that the 
ratio of the Higgs field inside the bubble to the temperature T cannot be 
smaller than one, in order that the sphaleron energy not be too small. 
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This condition was used in ITT], |T9l to impose a strong constraint on the 



Higgs mass in the minimal version of the electroweak theory, mn ^ 42 GeV. 
This, of course, aheady contradicts the present experimental limits mn ^ 57 
GeV [EU|. However, more careful consideration of various theoretical uncer- 



tainties indicated that the constraint might be somewhat weaker, permitting 
rriH up to 55 GeV, or possibly higher pT|. In any case, successful baryogene- 



sis almost certainly requires some extension of the standard model, possibly 
including more scalar fields. In multi-Higgs models |19, the limits are 



substantially weaker. Indeed, Anderson and Hall p2[ have noted that simply 



adding a scalar singlet to the model significantly weakens the constraint. 

Before one can discuss details of the process of baryogenesis, it is necessary 
to check that the results of our investigation of the phase transition are 
reliable. This is not a trivial issue even in the minimal electroweak theory. 
Indeed, as stressed in Refs. p, |^ , each new order of perturbation theory at 
finite temperature may bring a new factor of g'^T/m ~ gT / cj) for the theories 
with gauge boson masses m ~ gcf). This means that the results of the one 
loop calculations may become unreliable at ^ gT. A rather unexpected 
(and often overlooked) consequence of this observation is that we cannot 
even say in an absolutely reliable way that symmetry at high temperatures 
is completely restored; one may say only that the strength of the symmetry 
breaking is limited by the constraint < gT, m < g^T 0. In theories with 
g I this result is quite informative. However, gauge coupling constants 
in the electroweak theory are not much smaller than one. Therefore, the 
reliability of our results concerning the region (p ^ T deserves a detailed 
investigation. 

This issue became more urgent with the recent claim by Brahm and Hsu 
that higher order corrections lead to the appearance of a term in the effective 
potential ~ —g^(j)T^ [^. This term is linear in 0; it is very large at small 



0, and it removes the local minimum of V{(j), T) created by the cubic term 
~ ~g^(j)^T . As a result, the phase transition ceases to be first order. In 
fact, the phase transition ceases to be a phase transition, since the scalar 
field does not vanish at any temperature. At large temperature the scalar 
field remains smaller than gT, which means that this result is not in direct 



contradiction with our earlier conclusions. But still, the result of Ref. [g3 
seems somewhat surprising. 
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The importance of understanding these issues has been highhghted even 
more dramatically by the recent work of Shaposhnikov |^. He has also 



found a linear term in the effective potential, but with an opposite sign! 
Shaposhnikov concludes that the phase transition is much more strongly first 
order than expected. He argues that baryogenesis in the minimal electroweak 
model is possible despite the problems with the CP violation and obtains an 
improved constraint on the Higgs boson mass, rriH ^ 64 GeV, which is quite 
consistent with the experimental constraints. Thus, without a proper study 
of the higher order corrections to the effective potential, one may be unable 
to make any conclusions concerning the possibility of baryogenesis in the 
standard model. 



The authors of Refs. [psj |24| obtained linear terms by simply substitut- 
ing the effective masses found at one loop back into the one loop calculation 
of T). Such a procedure is generally reliable when calculating Green's 
functions, or tadpoles corresponding to ^(0, T). However, it leads to com- 
binatoric errors when calculating the free energy. Our investigation of this 
problem shows that if one is careful with counting of Feynman diagrams and 
with gauge invariance, neither positive nor negative linear terms ~ g^(f)T^ 
appear in the effective potential [p5l . 

Moreover, we will argue that, despite all uncertainties with higher order 
corrections, the expectation value of the scalar field actually disappears, 
= 0, at a temperature higher than some critical temperature To. (Note, 
that this would be impossible in the presence of linear terms of either sign.) 

However, high order corrections do lead to a definite and significant mod- 
ification of the one-loop results. Namely, they lead to a decrease of the cubic 
term g^cfr'T by a factor 2/3. This effect decreases the ratio 0/T at the point 
of the phase transition by approximately the same factor 2/3. This makes 
baryogenesis virtually impossible in the context of the minimal standard 
model with mn > 57 GeV. 

Assuming that one knows the shape of the effective potential at small 0, 
one should still work hard to determine the ratio 0/T at the point of the 
phase transition. One needs to know at what temperature the transition ac- 
tually occurs, and some details of how it occurs. At very high temperatures 
the effective potential of the Higgs field, V{(f),T), has a unique minimum 
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at the symmetric point (f) = 0. As the temperature is lowered, a second 
minimum appears. At a critical value Tc, this second minimum becomes de- 
generate with the first one. However, the phase transition actually occurs 
at a somewhat lower temperature, due to the formation of bubbles of true 
vacuum which grow and fill the universe. The usual way to study bubble 
formation is to use the euclidean approach to tunneling at a finite temper- 



ature [^]. One should find high-temperature solutions, which describe the 
so-called critical bubbles. Then one should calculate their action, which leads 
to an exponential suppression of the probability of bubble formation. Typi- 
cally, these calculations are rather complicated, and analytic results can only 
be obtained in a few cases. One of these is the thin wall approximation, 
which is valid (as in the case of transitions at zero temperature) if the differ- 
ence in depth of the two minima of V{(f), T) is much smaller than the height 
of the barrier between them. In this case the radius of the bubble at the 
moment of its formation is much larger than the size of the bubble wall, and 
the properties of the bubble can be obtained very easily. Recently Ander- 
son and Hall performed a thorough analytic study of the phase transition in 



the electroweak theory in one loop approximation ||22|. Their results for the 
one loop effective potential V{(f),T) completely agree with the results of our 
investigation |21[]. However, in their study of the bubble formation they as- 



sumed that the thin-wall approximation is applicable. As we will see, for the 
Higgs boson masses in the range of interest, ttih < this is not the case 

even if one takes into account modification of the cubic terms by higher order 
corrections. Fortunately, we were able to obtain a simple analytic expression 
which gives the value of the euclidean action for theories with an effective 
potentials of a rather general type, V^(0, T) = — bcp^ + ccp^. We hope that 
this result will be useful for a future investigation of bubble formation in a 
wide class of gauge theories with spontaneous symmetry breaking. 

On the other hand, validity of the standard assumption that the phase 
transition occurs due to formation of critical bubbles should be verified as 
well. Kolb and Gleiser [53] and, more recently, Tetradis P9 have argued 



that the phase transition may occur by a different mechanism, the formation 
of small (subcritical) bubbles. If this is the case, the transition is completed 
earlier and by a different mechanism than in the conventional picture. While 
this idea is very interesting, we will argue that it is only relevant in cases 
where the transition is very weakly first order and the euclidean action cor- 
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responding to critical bubbles is not much larger than one. This is not the 
case for the strongly first order phase transitions, where the relevant value 
of the euclidean action at the moment of the transition is 5 ~ 130 — 140. 

Determination of the baryon asymmetry produced at the phase transi- 
tion requires knowledge not only of how bubbles are produced, but also of 
how they evolve. Because the expansion of the universe is so slow at this 
time, a typical bubble grows to a macroscopic size before colliding with other 
bubbles. In the first scenarios proposed for the formation of the asymmetry, 
baryon number was produced in the bubble wall |jll| - |T^. This mecha- 
nism, at best, is not terribly efficient, because the baryon number violating 
processes turn off rapidly as the scalar field expectation value turns on. The 
resulting asymmetry is sensitive to the speed and thickness of the bubble.^ 
The most effective scenarios for electroweak baryogenesis have the baryons 



produced in front of the wall, in the symmetric phase ||T8[. In this pic- 
ture, scattering, for example, of top quarks from the bubble wall leads to an 
asymmetry in left vs. right-handed top quarks in a region near the wall. This 
asymmetry, resulting from an asymmetry between reflection and transmis- 
sion of different quark helicities at the wall, biases the rate of baryon number 
violation in the region in front of the wall; the resulting value of nfe/n^ can be 



as large as 10~ . However, the authors of Ref. [|I^ assumed that the wall was 
rather thin, with a thickness of order T~^. For thicker walls, the asymme- 
try goes rapidly to zero. This can easily be understood. In order to have an 
asymmetry in reflection coefficients, the top quarks must have enough energy 
to pass through the wall. For rrit ~ 120 GeV, this means typically the energy 
must be greater than about T/2. If the wall is very thick compared to this 
scale, the motion of the top quarks is to a good approximation semiclassical, 
and the reflection coefficient is exponentially suppressed. The analyses of 
other authors also exhibit sensitivity to the wall shape and velocity. 

Clearly, then, it is important to understand how the bubble propagates 
after its initial formation. In this paper we will consider some aspects of this 
problem. A complete description of the wall evolution is rather complicated. 
We will see, however, that in certain limits it is not too difficult to determine 
how the velocity and thickness of the wall depend on the underlying model 



""^The extent of this suppression and the precise form of this dependence is still a subject 
of debate. 
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parameters. In this analysis, it is crucial to recall that in addition to the vari- 
ous microscopic parameters, there is another parameter of great importance: 
the expansion rate of the universe. Imagine a world with arbitrarily small 
value of Newton's constant. In such a world, the phase transition occurs at 
a temperature arbitrarily close to and the pressure difference on the two 
sides of the wall is arbitrarily small. Thus, in this limit one expects that the 
wall will move arbitrarily slowly, and a systematic expansion of the relevant 
physical quantities in powers of the wall velocity should be possible. Since in 
the real world, the expansion rate of the universe is indeed quite small at Tc, 
it seems plausible that the real velocity of the wall should be small, and that 
the expansion in powers of velocity should be a fair approximation. Even 
within the framework of this small velocity approximation, we will content 
ourselves with a quite crude picture for the processes which damp the wall's 
motion. Within this framework, we will find that for plausible values of Higgs 
masses, the wall velocity is indeed small, varying from v ~ .01 to w ~ 0.3; at 
the same time the wall will typically have a thickness of order lO's of T~^'s. 

There have been a number of attempts in the past to compute the wall 
parameters, including efforts by some of the present authors. In Ref. 
a simple formula for the wall velocity was given, based on a semiclassical 
picture in which one species of particles gains a large mass M ^ T as 
it passes through the wall. Balancing the force on the wall due to these 
particles with the pressure difference between the two phases gives a relation 
of the form 

Ap 

where p is the pressure and p is the internal energy. We will see that an 
expression of this form holds only if the mean free paths of the particles are 
large compared to the thickness of the wall, and that this condition is not 
likely to hold for the electroweak transition. We will argue that the actual 
wall velocity is somewhat larger than would be expected from this formula. 



In Ref. |]T5[ , a different approach was adopted. Assuming the temperature 
and velocity are constant across the wall, momentum and energy conservation 
give an expression for the velocity of the form 
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As we will see in the present investigation, however, one cannot neglect the 
change in temperature and velocity of the gas as it passes the wall. For 
relativistic gases, one always obtains expressions linear in the velocity. 



More recently, Turok has argued that these types of analyses are in- 
correct. He suggests that reflection of particles from the wall does not slow 
the wall at all. Indeed, he argues that the only way in which the wall can 
dissipate energy is through rather complicated particle production processes, 
suppressed by several powers of coupling constant. As a result, even if the 
phase transition is very weakly first order, the bubble wall becomes ultrarel- 
ativistic. To buttress his case, Turok shows that if the gas is everywhere in 
thermal equilibrium (corresponding to the local value of the scalar field) then 
the force on the wall is independent of the velocity. This is indeed correct. 
However, all of the effects we find arise because of small, velocity-dependent 
departures from equilibrium.^ 

To a large extent, our analysis will be an extension of the semiclassical 



reasoning of Ref. ||2^. We begin by pointing out that there is an important 
assumption made more or less explicitly in all of these analyses: that by 
the time the wall has grown to a macroscopic size, it has achieved a steady 
state. If this is the case, the problem can be analyzed in the rest frame of 
the wall, where the scalar field and the particle distributions are independent 
of time. We describe a simple situation (due to L. Susskind) where this is 
not the case, and which we refer to as the 'snowplow.' In this situation, 
there is a steady pileup of particles near the wall. While we do not expect 
that precisely this snowplow phenomenon occurs in the cases of interest, 
it makes clear that there are additional, potentially important effects which 
must be taken into account and which are left out of existing treatments. We 
then consider there three limiting situations. In the first of these, a typical 
particle does not scatter off of other particles of the gas as it passes through 
the wall. This 'thin wall' case is that for which (a modification of) eq. (j^) 
is valid, and tends to give a very small value for the wall velocity. In the 
second case, the wall is thick compared both to typical mean free paths for 
elastic scattering and for scatterings which change the numbers of different 



^Turok and McLerran have recently informed us that they have reconsidered the ar- 
gument of ||30[| and have argued that there are additional sources of damping, including 
scattering of particles from the wall. 
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particle types. In this case, the gas is nearly in equilibrium everywhere. The 
velocity is larger than in the thin wall case by a factor of order y^, where 
^ is some typical mean free path and 5 is the wall thickness. In the third 
situation, the wall is thick compared to typical mean free paths for elastic 
scattering, but not for scatterings which change particle number. Thus, one 
has some approximate kinetic equilibrium locally, but particle numbers are 
not equilibrated in the wall. This is the situation which appears to have 
greatest relevance to the electroweak transition. Here, phenomena similar 
to the snowplow effect occur, and there is an enhancement of the density of 
tops, W's and Z's in the wall. This tends to reduce the wall velocity, giving 
a result intermediate between the thin and thick wall cases. 

The plan of the paper is the following. In Section 2 we will describe the 
phase transition in the electroweak theory in one loop approximation. In 
Section 3 we will consider the theory of bubble formation during the phase 
transition. Section 4 will contain a discussion of role of the higher order 
corrections. In Section 5 we will deal with the issue of subcritical bubbles. 
Finally, in Section 6, the bubble wall propagation will be considered. Details 
of relevant calculations will be contained in the Appendix. 



2 The Phase Transition 



Let us consider the form of the effective potential at finite temperature. 
Contributions of particles of a mass m to 1^(0, T) are proportional to T^, 
im? T and ln(m/T). We will assume that the Higgs boson mass is smaller 
than the masses of W and Z bosons and the top quark, mn < t^w^ ^t- 
Therefore we will neglect the Higgs boson contribution to V{(f),T). 

The zero temperature potential, taking into account one-loop corrections, 
is given by |^ 

Vo = -^<P' + j^' + 2Bvy - + B<p' In(^) . (3) 

Here 

i? = ^(2m^^ + m|-4mn, (4) 
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Vo = 246 GeV is the value of the scalar field at the minimum of Vq, A = 
/i^/f^, rriff = 2/i^. Note that these relations between X,fi,Vo and the Higgs 
boson mass mn, which are true at the classical level, are satisfied even with 
an account taken of the one-loop corrections. This is an advantage of the 
normalization conditions used in |^ . An expression used in |2^ is equivalent 
to this expression up to an obvious change of variables. 

At a finite temperature, one should add to this expression the term 

Vt = ^{6i-{yw) + ^i-iyz) - Gi+{yt)) , (5) 

where i/i = MiCp/voT, and 

Uy) = ± r dx ln(l T e~^^^) . (6) 
Jo 

The results of our work are based on numerical calculation of these integrals, 
without making any specific approximations |21]. However, in the large tem- 
perature limit it is sufficient to use an approximate expression for V{(j), T) 



Here 



T) = D{T^ - T^)(j)^ - ET(f + . (7) 

D = -^(2m2^ + m| + 2m?) , (8) 



E = ^(2m3^ + m|)~10-2, (9) 

T'o = ^(/i^ - ^Bvl) = ^{mj, - Wvl) , (10) 

At = A ^f2mti.ln^ + miln^-4m^ln^^ , (11) 

where In = 2 ln47r — 27 ~ 3.91, InOi? = 2 Invr — 27 ~ 1.14. 

To avoid misunderstandings, we should note again that due to our choice 
of more convenient renormalization conditions, the form of some of our equa- 
tions is slightly different from the form of expressions used in [Q. In partic- 
ular, instead of large coefficients cb and cp in the equation for At in [^], we 
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have smaller constants as and ap. However, all physical results obtained by 
our equations coincide with the one-loop results of Refs. . 

The behavior of V{(f),T) is reviewed in Refs. 0. It will be useful for 
our future discussion to identify several 'critical points' in the evolution of 
Vi<jy,T). 

At very high temperatures the only minimum of V{(j),T) is at </> = 0. A 
second minimum appears at T = Ti, where 

t2 ^ ±o /12^ 

^ 1 - dEysXr.D ■ ^ ^ 



The value of the field (f) in this minimum at T = Ti is equal to 

3ETi 



2A 



Ti 



(13) 



The values of V{(f), T) in the two minima become equal to each other at the 
temperature Tc, where 

'Ji2 

T2 = . (14) 

At that moment the field in the second minimum becomes equal to 

^0 = . (15) 

The minimum of V{(j),T) at = disappears at the temperature To, when 
the field cj) in the second minimum becomes equal to 

= ^ . (16) 

The results of a numerical investigation of V{(j), T) for a particular case, 
irtH = 50 GeV and mt = 120 GeV, are shown in Fig. 1. 



3 Bubble Formation 



In the previous section we noted that the two minima of V{(f), T) become of 
the same depth at the temperature eq. (|14|). However, tunneling with 
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formation of bubbles of the field cj) corresponding to the second minimum 
starts somewhat later, and it goes sufficiently fast to fill the whole universe 
with the bubbles of the new phase only at some lower temperature T when 
the corresponding euclidean action suppressing the tunneling becomes less 
than 130 - 140 |]13|, |21|, Some small uncertainty in this number is related 



to the speed with which bubble walls move after being formed (see the next 
section) and to the exact value of the critical temperature, which is very 
sensitive to the top quark mass and to the value of the cubic term. In this 



paper (see also |^) we performed a numerical study of the probability of 
tunneling. Before reporting our results, we will remind the reader of some 
basic concepts of the theory of tunneling at a finite temperature. 



In the euclidean approach to tunneling (at zero temperature) the 
probability of bubble formation in quantum field theory is proportional to 
exp(— 5*4), where 5*4 is the four-dimensional Euclidean action correspond- 
ing to the tunneling trajectory. In other words, 5*4 is the instanton action, 
where the instanton is the solution of the euclidean field equations describ- 
ing tunneling. A generalization of this method for tunneling at a very high 



temperature gives the probability of tunneling per unit time per unit 
volume 

P~A(T)-exp(-^). (17) 

Here A{T) is some subexponential factor roughly of order T^; 5*3 is a three- 
dimensional instanton action. It has the same meaning (and value) as the 
fiuctuation of the free energy F = V^(0(x),T) which is necessary for bubble 
formation. To find S3, one should first find an 0(3)-symmetric solution, (j){r), 
of the equation 

^ + = (18) 

with the boundary conditions 0(r = 00) = and dcf) / dr\r=o = 0. Here 

r = \fx^; the Xi are the euclidean coordinates, i = 1,2,3. Then one should 
calculate the corresponding action 

^3 = 47r 1^°° dr[^ (^^^ + V(0(r), T)] . (19) 
Usually it is impossible to find an exact solution of eq. ( |l8l) and to 
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calculate 5*3 without the help of a computer. A few exceptions to this rule 
are given in Refs. 0, |2^. One of these exceptional cases is reahzed if the 
effective potential has two almost degenerate minima, such that the difference 
e between the values of V{(f),T) at these minima is much smaller than the 
energy barrier between them. In such the thickness of the bubble wall 

at the moment of its formation is much smaller than the radius of the bubble, 
and the action 5*3 can be calculated exactly as a function of the bubble radius 
r, the energy difference AV and the bubble wall surface energy 5*1: 

All 

S; = -—r''AV + A'Kr^Si, (20) 

o 

where 

5i = d<p^2V{^,T) . (21) 

The radius of the critical bubble r can be found by finding an extremum of 
5*3 (r). However, one must be very careful when using these results. Indeed, 
as can be easily checked, this extremum is not a minimum of the action, it 
is a maximum. (This just corresponds to the fact that critical bubbles are 
unstable and either expand or contract). Similarly, the action correspond- 
ing to the true solution of eq. ([T8| ) will be higher than the action of any 
approximate solution. As a result, one can strongly overestimate the tun- 
neling probability by calculating it outside the limit of validity of the thin 



wall approximation. For example, in Ref. |2^ the phase transition in the 
electroweak theory with rriH ~ 50 GeV was studied and it was found that 
it happens very soon after the temperature approaches Tc, occurring due to 



formation of bubbles with thin walls. According to this happens at 
e = 1/6, where e = . However, the authors did not check the validity 

^ c ^ o 

of the thin wall approximation in this case. Whereas our one-loop results 
for the effective potential V{(f), T) are in complete agreement with the results 



of Ref. ||22[, our conclusion concerning the bubble formation is somewhat 
different. As we have already mentioned, the phase transition in the elec- 
troweak theory is completed when the ratio S3/T becomes about 130 - 140. 
Our calculations show that this happens at e ~ 1/4. In Fig. 1 we plot the 
shape of the effective potential and in Fig. 2 the shape of the solution of 
eq. ([Tsl ) corresponding to S3/T ~ 140 for rriH = 50 GeV and rrit = 120 
GeV. (The results for tunneling and for the ratio 0/T prove to be not very 
sensitive to the mass of the top quark in the interval 100 GeV ^ rrit ^ 150 
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GeV.) The effective potential V{(j),T) at e = 1/6 loolcs very similar, but the 
value of S3/T in this case is about 300, which is 3 times larger than the 



result obtained in |^2[ in the thin wall approximation. It is clear from these 
figures and the results of numerical calculations of S^/T that the thin wall 
approximation is far from being applicable for the investigation of the phase 
transition in the electroweak theory, unless the phase transition is weakly 
first order. However, the last case is not particularly interesting from the 
point of view of baryon asymmetry generation. 

We must note, that the numerical results obtained above are modified 
when one takes into account higher order corrections to the effective po- 
tential. As we will show in the next section, these corrections change the 
numerical value of the coefficient E in the cubic term in (|^. The final nu- 
merical results of our study of the probability of tunneling will be contained, 
therefore, in the next section. Here we just wanted to show the difference be- 
tween the results of the numerical investigation of tunneling and the results 
obtained in the thin wall approximation. This difference remains large after 
the modification of the coefficient E. 

We would now like to obtain an analytic estimate of the probability of 
tunneling in the electroweak theory, which can be used for any particular 
numerical values of constants D, E and At- As shown in Ref. |2^, eq. (0) 



in most interesting cases approximates V^(0, T) with an accuracy of a few 
percent. This by itself does not help very much if one must study tunneling 
anew for each new set of the constants. However, it proves possible to reduce 
this study to the calculation of one function f{a), where a is some ratio of 
constants D, E and A^- In what follows we will calculate this function for a 
wide range of values of a. This will make it possible to investigate tunneling 
in the electroweak theory without any further use of computers. 

First of all, let us represent the effective Lagrangian T) near the 
point of the phase transition in the following form: 

m T) = lid.^r - ^^0^ + w - . (22) 

Here M'^{T) = 2D{T^ — T^) is the effective mass squared of the field near 
the point = 0, Af is the value of the effective coupling constant At near 
the point of the phase transition (i.e. at T ~ Tt, where Tt is the temperature 
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at the moment of tunneling). With a very good accuracy, the constants 
Xt, Xti, are equal to each other. 



as: 



Defining = x = X/M, the effective Lagrangian can be written 



L{^,T) 



2^ ^ ^ 2 2 8 



where 



a 



XqM^ 
2E2T2 



The overall factor -^^^ does not affect the Lagrange equation 

1 = $ $ + -«$ . 

dR^ RdR 2 2 



(23) 
(24) 

(25) 



Solving this equation and integrating over d^X = M ^d^x gives the following 
expression for the corresponding action: 



T 



4.85 

^2 r3 



X /(a) . 



(26) 



The function /(a) is shown in Fig. 3. It is equal |2^ to 1 at a = 0, and 
blows up when a approaches 1. In the interval from to 0.95 this function, 
with an accuracy about 1%, is given by the following simple expression: 



/(«) 



2.4 

\ — a 



0.26 



ay 



(27) 



In the vicinity of the critical temperature Tq, i.e. at AT = T — T^ <ti To, 
the action (p6D can be written in the following form: 

^3 38.8/^3/2 /AT\3/2 2XoDAT^ 

Using these results, one can easily get analytical expressions for the tunneling 
probability in a wide class of theories with spontaneous symmetry breaking, 
including GUTs and the minimal electroweak theory. 
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4 Infrared Problems and Reliability of the 
Perturbation Expansion 



In our previous discussion, we have considered only the one loop corrections 
to the effective potential. In this section we discuss the role of higher order 
corrections. 

Early investigations of the electroweak phase transition did not take 
into account corrections due to strong interactions, since at that time most 
physicists did not expect that top quarks would be heavier than W and Z 
bosons, and thus their contributions were not expected to be terribly impor- 
tant. Given what we now know about the top quark mass, it appears that 
top quarks give the largest contribution to the parameters D and in eqs. 
(H), (pTI) . Thus, it is natural to ask whether strong interaction corrections 
are likely to be important. 

Our preliminary investigation of this question indicates that this is not 
the case. For example, one of the most important effects would be a change 
of the Fermi distribution at a finite temperature due to the modification of 
the quark mass by interactions with gluons. According to quarks at a 
high temperature acquire a correction to their effective mass squared, 

AmfiT) = ^T' (29) 
6 

with Qs the strong coupling constant. ^ This gives Am^(T) ~ 0.2 at the 
temperature of the electroweak phase transition, T ~ 10^ GeV. A simi- 
lar contribution to the boson mass could lead to important effects (see be- 
low). However, due to Fermi statistics, fermion propagators contain terms 
[(2n + 1) vrT]^ 0. As a result, all thermodynamic quantities rather weakly 
depend on the effective mass squared of the fermions, until this mass becomes 
comparable with vrT. 

One should note, of course, that Am^(T) ~ 0.2 is not a true correc- 
tion to the top quark mass squared; rather it is a square of the mass gap in 
the spectrum of fermionic excitations. Therefore a more detailed analysis of 

■^We use Qfs ~ 0.1 for our estimates. 
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the higher-loop diagrams involving strong interactions is desirable. Never- 
theless, our estimate suggests that the strong interaction effects are actually 
insignificant for the description of the phase transition. They just lead to a 
small modification of the critical temperature. More importantly, they do 
not change cubic terms ~ g^(f)^T in the effective potential and do not induce 
linear terms ~ g^(f)T^. 

However, higher order corrections in weak interactions at T > Tc may 
be very important. It is well known that, in field theories of massless parti- 
cles, perturbation theory at finite temperature is subject to severe infrared 
divergence problems. For small values of the scalar field, the gauge bosons 
(and near the phase transition, the Higgs boson) are nearly massless; as a 
result, as was pointed out in the early work on this subject ^ |], one cannot 
reliably compute the effective potential for very small 0. In this section, we 
attempt to determine whether the standard one loop calculation of V is in- 
deed reliable in the range of relevant to our analysis. One might worry, for 
example, that since the term —ET(f)^ in the effective potential, which leads 
to the first order phase transition, is important only for rather small (f), there 
might be large corrections changing the order of the phase transition. We 
will show that indeed the coefficient of the (p^ term is altered to 2/3 of its 
one loop value. This renders the phase transition, for a given value of the 
couplings, less first order and can have significant effects on baryogenesis. 
On the other hand, we will argue that perturbation theory is not in terribly 
bad shape, and that one can determine the nature of the phase transition 
with some confidence from low order calculations. 

Recently, in a very interesting paper, Brahm and Hsu reach the opposite 
conclusion [^. These authors find that at small 0, higher order correc- 
tions to the scalar field contribution to the effective potential contain a large 
negative linear term —g^(j)T^, which eliminates any trace of a first order 
transition. They argue, moreover, that their calculation is reliable, i.e. that 
all other higher order corrections are under control and do not modify their 
conclusion. 

On the other hand, Shaposhnikov considers higher order corrections to 
the vector particle contribution to V{(f), T) and finds a large positive term 
+ g^(j)T^. He concludes that the phase transition is strongly first order 
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(0/T > 1) even for mn ~ 64 GeV |4| 



We will show that neither positive nor negative linear terms appear in 
the expression for V{(j), T) if one studies higher order corrections paying par- 
ticular attention to the correct counting of Feynman diagrams. We 



will employ two separate approaches to this problem: a straightforward enu- 
meration of Feynman diagrams, and an effective action analysis valid for a 
discussion of infrared effects. 

We will consider here for simplicity the contribution of the scalar particles 
and the W bosons only; adding the contribution of Z bosons is trivial. As 
we have already noted, for questions of infrared behavior, fermions may be 
ignored. Coulomb gauge, V ■ W = 0, is particularly convenient for the 
analysis, though the problem can be analyzed in other gauges as well. In 
this gauge, the vector field propagator D^,^ after symmetry breaking (and 
after a proper diagonalization) splits into two pieces, the Coulomb piece, Dqq, 
and the transverse piece, Dij. For non-zero values of the discrete frequency, 
ujn = 2TTnT, the Coulomb piece mixes with the 'Goldstone' boson. However, 
for the infrared problems which concern us here, we are only interested in 
the propagators at zero frequency. For these there is no mixing. One has 

Doo{uJ = 0,k)= \ (30) 



and 



A,(^ = 0, ^) = \ R,{k) , (31) 

where Pij = 5ij — The mass of the vector field W at the classical level is 
given by mw = QVo/'^- Propagators of the Higgs field and of the 'Goldstone' 
field X ill this gauge are given by 

dS) = rrr-^ ' (32) 

D^ik) = 1 . (33) 
k^ 

Let us review several ways of obtaining the standard one-loop expression 
for the cubic term in the effective potential, eq. (|^). The most straight- 
forward is to carefully expand eq. (|^) for the effective potential in yw = 
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^Ft ~ 2T- Indeed, the contribution of ly-bosons to the effective potential at 
T > mw{4>) is given by 



VW(0,T) = 




Here the expression in brackets coincides with the contribution of a scalar 
field with mass mw', the factor 2 appears since there are two VT-bosons with 
opposite charges, while the factor 3, which will be particularly important in 
what follows, corresponds to the two transverse and one longitudinal degrees 
of freedom with mass m\Y- 

Alternatively, we can obtain the cubic term by looking directly at the one- 
loop Feynman diagrams. For this purpose, it is only necessary to examine 
the zero frequency contributions. Certain diagrams containing four external 
lines of the classical scalar field naively give a contribution proportional to 
g^(j)'^] the cubic term arises because the zero frequency integrals diverge for 
small mass as T /m\^ ~ T / gcj). 

Consider, in particular, the zero frequency part of the expression for the 
one loop free energy in momentum space. It is simplest to compute the 
tadpole diagrams for dV/d(f), as indicated in Fig. 4, and afterwards integrate 
with respect to (p. The transverse gauge bosons give a contribution 

# 2 J (2i)3 + rnl, 8,r V »- • ^ ' 

where, by keeping only the zero frequency mode, we have dropped terms 
which are analytic in The Coulomb lines give half the result of eq. 

(P^l). Integration of the total vector field contribution correctly represents 
the cubic term in (|3^ . 

A complete gauge boson contribution to the tadpole, including the non- 

^ The integral has been defined by dimensional regularization; ultraviolet divergences 
here are absorbed into the usual zero temperature renormalizations. 
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zero frequency modes, is @] 



48 ^ TT ^ 48 ^ 277 ^ 

(36) 

One can easily check that integration of this expression with respect to 



gives eq. (|3^). 



With these techniques, we are in a good position to study higher order 
corrections to the potentiaL The authors of Refs. |2^, |Q found a hnear 



contribution to the potential by substituting the mass found at one loop back 
into the one loop calculation. The effective masses-squared of both scalar 
particles and of the Coulomb field contain terms of the form ~ g'^T(j), which, 
upon substitution in (34), give linear terms. But this procedure is not always 
correct. It is well known that the sum of the geometric progression, which 
appears after the insertion of an arbitrary number of polarization operators 
n((/), T) into the propagator (Z^^ + m^)"^, simply gives [k"^ + m'^ + Il{(j),T))~^. 
Therefore one can actually use propagators {k'^ + m'^+Il{T))^^, which contain 
the effective mass-squared + 11(0, T) instead of m^. However, this trick 
with the geometric progression does not work for the closed loop diagram 
for the effective potential, which contains ln(fc2 -|- m^). A naive substitution 
of the effective mass squared + 11(0, T) instead of into ln{k^ + m^) 
corresponds to a wrong counting of higher order corrections. 

A simple way to take into account high temperature corrections to masses 
of vector and scalar particles without any problems with combinatorics is to 
compute tadpole diagrams for ^; these are then trivially integrated to give 
the potential. One can easily check by this method that no linear terms 
appear in the expression for 1^(0, T). Indeed, at a given temperature and 
effective mass, the tadpoles are linear in (see e.g. equation (^)). To take 
into account the mass renormalization in the tadpoles, one should substitute 
the effective mass squared -|- 11(0, T) into the one-loop expression for 
the tadpole contribution; as we explained above (see also (0)), this is a 
correct and unambiguous procedure for tadpoles. Since + 11(0, T) is not 
singular in the limit — > 0, the tadpole ( PB] ) in this limit remains linear in 
0. Therefore its integration with respect to 0, which gives the correction to 
the effective potential, is quadratic in 0, i.e. it does not contain any linear 
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terms. 



Even though there are no hnear terms ~ g^(f)T^, higher order corrections 
do have a dramatic effect on the phase transition, which has apparently not 
been noted before. This effect is a modification of the cubic term. 

As we have shown above, the cubic term appears due to the contribution 
of zero modes, u;„ = 27mT = 0. This makes it particularly easy to study 
its modification by high order effects. Indeed, it is well known that the 
Coulomb field at zero frequency acquires the Debye 'mass', m^, = noo(co'n = 
0, fc — > 0) ~ g'^T'^- This leads to an important modification of the Coulomb 
propagator (PU]): 

^ = 2^ 2u^ ■ (^^^ 

For the values of (f) of interest to us, mj^ ^ m^(0). Thus, repeating the 
calculation of the cubic term, the Coulomb contribution disappears. How- 
ever, the transverse contribution, which is two times larger than the Coulomb 
one, is unaffected at this order, due to the vanishing of the 'magnetic mass' 



^The absence of a linear term in </> can be understood in an effective Lagrangian ap- 
proach as well. Such an approach automatically gives the correct combinatorics. Since, 
at weak coupling, we are interested in energy scales much less than ttT, one can first 
integrate out all modes of the various fields with non-zero Matsubara frequency, thereby 
obtaining a three-dimensional effective action for the light fields. At one loop, one has 
the usual mass corrections, order g^T^ for the Coulomb and scalar lines, whereas the 
quadratic term for the transverse gauge bosons goes to zero with k and (f>. These 'mass 
corrections' are analytic in in particular, they do not contain linear terms in ip. By 
assumption, the Higgs field is light at the phase transition, but one can again integrate 
out the massive Coulomb field. It is still true (as a consequence of gauge invariance) that 
the quadratic term for the transverse gauge bosons vanishes at zero momentum. Having 
done this procedure, we compute the effective potential for (f) by calculating the determi- 
nant of Gaussian fiuctuations in the low energy theory, i.e. the three dimensional theory 
containing only transverse gauge bosons and scalars. The linear terms in (p used in Ref. 

come from loops of light fields (lu = 0), and so in the present language we need to 
check that there is no non-analytic behavior in higher loop graphs in the effective theory. 
In fact, although individual diagrams may be singular, the sum of all two loop graphs is 
non-singular, essentially due to gauge invariance, and so the effective potential contains 
no linear term in (/), to this level of approximation. The cubic term is correctly reproduced 
by this process (see below). 
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[|, As a result, the cubic term does not disappear, but it is diminished 
by a factoif]of 2/3: [] 

E = -^^{2m'^ + ml) . (38) 

This small correction proves to be very significant. Indeed, eqs. ([T5|), (|1^) 
show that the ratio of the scalar field to the temperature at the moment 
of the phase transition is proportional to E, i.e. to the cubic term. Actually, 
the dependence is even slightly stronger, since for smaller E the tunneling 
occurs earlier. Results of a complete numerical investigation of the ratio 
(j)/T at the moment of the phase transition, as a function of the Higgs boson 
mass are shown in Fig. 5, for the top quark mass rrit = 120 GeV. We have 
found that the ratio 0/T is not very sensitive to the mass of the top quark, 
in the interval 100 GeV ^ rrit ^ 150 GeV, and it decreases for rrit outside 
this interval. Even before the reduction of the cubic term was taken into 
account, the ratio 0/T for ttih ^ 57 GeV was shghtly less than the critical 
value 0/T ^ 1. The decrease of this quantity by a factor of 2/3 makes it 
absolutely impossible to preserve the baryon asymmetry generated during 
the phase transition in the minimal model of electroweak interactions with 
rriH ^ 57 GeV. 

Is this the end of the story? The effective coupling constant of interac- 
tions between W bosons and Higgs particles is g/2. In this general 
investigation of the infrared problem in the non-Abelian gauge theories at a 
finite temperature suggests that the results which we obtained are reliable 
for ^ I T ~ T/3 Thus, a more detailed investigation is needed 



to study behavior of the theories with mn ^ 10^ GeV near the critical tem- 



^We should point out that if the magnetic mass of the transverse gauge bosons were 
to be numerically rather large, then the factor E would be reduced further, leading to a 
further weakening of the order of the phase transition. 

^ After we obtained this result, we received a paper by Carrington |Q where the renor- 
malization of the cubic term was also considered. A similar investigation was performed 
by Shaposhnikov as well |Q. These authors did not point out that they obtained the 
renormalization of the cubic term by the factor of 2/3, but after some algebra one can 
identify those terms in their expressions which are equivalent to ours. However, for differ- 
ent reasons, the final numerical results for the ratio (/)/r obtained by both Shaposhnikov 
and Carrington differ from our result considerably, being approximately two times larger 
than our result shown in Fig. 5. 
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perature, since the scalar field, which appears at the moment of the phase 
transition in these theories, is very small (see Fig. 5). However, we expect 
that our results are reliable for strongly first order phase transitions with 
(p ^ T, which is quite sufficient to study (or to rule out) baryogenesis in the 
electroweak theory. 

Finally, we would like to address a fundamental question: since the the- 
ory for (f) <^ gT is infrared divergent, can we definitely establish that the 
symmetry is restored at high temperature, or is it possible that always 
has some small, non-zero value? To address this question, we can work at 
T ^ To. In this case, the scalar field is massive, and scalar loops are not 
singular in the infrared. Potential infrared problems arise only from gauge 
boson loops. For these, the situation is similar to that in high tempera- 
ture QCD |3^. The standard assumption about QCD is that the infrared 
divergences are cut off at a scale rrimag ~ g^T (the detailed mechanism of 
the infrared cutoff will not be important to us). The free energy, Q, is non- 
singular through order g^T^. At order g^T"^, there is a logarithmic infrared 
divergence; Q ~ T^ln{mmag/T) ~ g^T^lng^. Higher order corrections go 
as g^T^{g'^T /mmag)^^ i-^- they are all of the same order. It is usually as- 
sumed, then, that the free energy can be computed through order g^T'^lng'^. 
A similar investigation suggests that the effective mass can be calculated 
with an accuracy g'^T^lng^. 

What are the consequences of these assumptions for the electroweak the- 
ory? First, the existence of an infrared cutoff of order g^T means that the 
potential is analytic in for small 0. This implies, in particular, that 
there are no linear terms in 0. Consider, then, the calculation of the 0^ 
term, which determines the value of the critical temperature. At lowest 
order, one has the standard result, DT'^cf? ~ g'^T'^cfP'. At two loop order, 
there is a correction proportional to g^T"^ (j)"^ In g"^; higher order corrections 
all go as g'^T'^cfP' . Thus, as in QCD, we can say that the mass term can be 
calculated to order g'^T'^lng'^. This is a small correction, and the one- loop 
calculation is reliable: for small 0, the curvature of the potential is the sum 
of the (zero-temperature) negative term — ^(f)"^ + 2BVo(p'^ , and a positive term 
DT'^cj)^ + 0{g^\ng'^)T'^(lP' which grows with temperature. Thus, we have a 
phase transition, and the higher loop effects can only lead to corrections 
~ g^ In g^ to the value of the critical temperature obtained in the one loop 
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approximation. 

Note that this discussion is vahd for any value of the Higgs mass. If these 
arguments are correct, then we expect the situation with the phase transitions 
in the non-Abehan gauge theories to be the same as in the standard case: 
infrared problems may prevent a simple description of the phase transition 
in a small vicinity of the critical point (unless the phase transition is strongly 
first order), but everywhere outside this region, the symmetry behavior of 
gauge theories can be described in a rehable way. We hope to return to a 
discussion of this interesting question in a separate publication. 



5 Subcritical Bubbles 



Despite our semi-optimistic conclusions concerning the infrared problem, it 
is still desirable to check that the whole picture of the behavior of the scalar 
field described above is (at least) self-consistent. This means that if the 
effective potential is actually given by eqs. (|^), (|]), (p^, (|TTD , (0), then 



our subsequent description of the phase transition and the bubble formation 
is correct. Indeed, one would expect that the theory of bubble formation 
is reliable, since the corresponding action for tunneling S^/T is very large, 
S3/T 130 — 140. However, recently even the validity of this basic assump- 
tion has been questioned. Gleiser and Kolb and Tetradis have argued 



that in many cases phase transitions occur not due to bubbles of a critical 
size, which we studied in section 3, but due to smaller, subcritical bubbles. 
We believe that these authors raise a real issue. However, we will now argue 
that this problem only arises if the phase transition is extremely weakly first 
order. 



The basic difference between the analysis of Ref. [p8| , |29| and the more 
conventional one is their assumption that at the time of the phase transition 
there is a comparable probability to find different parts of the universe in 



either of the two minima of V{(j),T). The main argument of Ref. ^ 
is that if the dispersion of thermal fluctuations of the scalar field < cfP' > 

is comparable with the distance between the two minima of V^(0, T), 
then the field d) "does not know" which minimum is true and which is false. 
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Therefore it spends comparable time in each of them. According to a 
kind of equihbrium between the domains of the two types is achieved due to 
subcritical bubbles with small action S^/T if many such bubbles may appear 
within a horizon of a radius H~^. 

In order to investigate this question in a more detailed way, let us re- 
examine our own assumptions concerning the distribution of the scalar field 
(f) prior to the moment at which the temperature drops down to Ti, when 
the second minimum of T) appears. According to (|13|), the value of the 
scalar field in the second minimum at the moment when it is formed is equal 
to 01 = For rriH ~ 60 GeV (and taking into account the coefficient 2/3 
in the cubic term) one obtains 0i ~ 0.4 T. Thermal fluctuations of the field 
have the dispersion squared < cfP' >= T^/12. (Note an important factor 
1/12, which was absent in the estimate made in [Q.) This gives dispersion 
of thermal fiuctuations ■\/< 0^ > ~ 0.3 T, which is not much smaller than 



However, as the authors of p8| emphasized in their previous work |^ 
(see also |2^), the total dispersion < 0^ >~ T^/12 is not an adequate 



quantity to consider since we are not really interested in infinitesimally small 
domains containing different values of fluctuating fleld 0. They argue that 
the proper measure of thermal fluctuations is the contribution to < 0^ > 
from fluctuations of the size of the correlation length ^(T) ~ M~^(T). This 
leads to an estimate < 0^ >~ T M(T), which also may be quite large [P^]. 
Here again one should be very careful to use the proper coefficients in the 
estimate. One needs to understand also why this estimate could be relevant. 



In order to make the arguments of Ref. [^, 29 1 more quantitative and 
to outline the domain of their validity, it is helpful to review the stochastic 
approach to tunneling (see |]3^ and references therein). This approach is 



not as precise as the euclidean approach (in theories where the euclidean 
approach is applicable). However, it is much simpler and more intuitive, 
and it may help us to look from a different point of view on the results we 
obtained in the previous section and on the approach suggested in 



The main idea of the stochastic approach can be illustrated by an example 
of tunneling with bubble formation from the point = in the theory (^) 



27 



with the effective potential 

V^(0,T) = ^^0^-ET0^ + ^0^ (39) 

For simphcity, we will study here the limiting case Aq — > 0. 

At the moment of its formation, the bubble wall does not move. In the 
limit of small bubble velocity, the equation of motion of the field at finite 
temperature is simply, 

= d^(j)/dr'^ + {2/r)d(j)/dr - 1/'(0) . (40) 

The bubble starts growing if > 0, which requires that 

\d^(t)/dr^ + {2/r)d<i)/dr\ < -^'(0) . (41) 

A bubble of a classical field is formed only if it contains a sufficiently big 
field 0. It should be over the barrier, so that dV/dcf) < 0, and the effective 
potential there should be negative since otherwise formation of a bubble will 
be energetically unfavorable. The last condition means that the field inside 
the critical bubble should be somewhat larger than 0*, where \^(0*,T) = 
V{<d,T). In the theory (||) with 0, one has 0* = M^/2ET. As 

a simplest (but educated) guess, let us take ~ 20^, = M^/ET. Another 
important condition is that the size of the bubble should be sufficiently large. 
If the size of the bubble is too small, the gradient terms are bigger than the 
term |V^'(0)|, and the field inside the bubble does not grow. Typically, the 
second term in ( PD somewhat compensates the first one. To make a very 



rough estimate, one may write the condition (O) in the form 



~ < ^fcL. - r'lV'm ~ (42) 

Let us estimate the probability of an event in which thermal fluctuations 
with T ^ M build up a configuration of the field satisfying this condition. 
The dispersion of thermal fluctuations of the field with k < k^ax is given 
by 

2 1 f'^™^ Pdk 

< >k<k 



VF+M^ f exp ^^^^^ 



T f^rnax k^ dk 

(43) 



27r2 Jo fc2 + M2 
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Note that the main contribution to the integral is given by k'^ ~ k'^ax ~ 4M2. 
This means that one can get a reasonably good estimate of < 0^ >k<kmax 
by omitting in the integrand. This also means that this estimate will be 
good enough even though the effective mass of the scalar field M'^[(j)) = V"{(j)) 
changes between = and (p. The result we get is 

T /-fcrnax Tk C'^TM 
< r >k<k^a.- TT^ dk= = — . 44 

Here C = 0(1) is a coefficient reflecting the uncertainty in the determination 
of kmax and estimating the integral. 

Thus, we have a rough estimate of the dispersion of perturbations which 



may sum up to produce a field which satisfies the condition (|42D . We can 
use it to evaluate the probability that these fluctuations build up a bubble 
of the field of a radius r > k^l^^. This can be done with the help of the 
Gaussian distribution]^ 

(45) 

Note that the factor in the exponent in ( P5| ) to within a factor of = 1.02 
coincides with the exact result for the tunneling probability in this theory 
obtained by the euclidean approach (see eq. (p^)): 

4 85M^ 

P - «-P(-^|^) • (46) 

Taking into account the very rough method we used to calculate the disper- 
sion of the perturbations responsible for tunneling, the coincidence is rather 
impressive. 



As was shown in ||3^ , most of the results obtained in the tunneling theory 
by euclidean methods can easily be reproduced (with an accuracy of the 
coefficient = 0(1) in the exponent) by this simple method. 

Now let us return to the issue of subcritical bubbles. As we have seen, 
dispersion of the long- wave perturbations of the scalar field, < 0^ >i.^k^^^c:^ 

^The probability distribution is approximately Gaussian even though the effective po- 
tential is not purely quadratic. The reason is that we were able to neglect the curvature 
of the effective potential = V" while calculating < (jp >fc<fc„„^- 
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is quite relevant to the theory of tunnehng. Its calculation provides 
a simple and intuitive way to get the same results as we obtained earlier by 



the euclidean approach [^. To get a good estimate of the probability of 



formation of a critical bubble in our simple model one should calculate this 
dispersion for k^ax ~ 2M(T), which gives < 0^ TM/vr^. Note, 

that this estimate is much smaller than the naive estimate < cfP >^ TM. 

The crucial test of our basic assumptions is a comparison of this dispersion 
and the value of the field at the moment T = Ti, when the minimum at 
= 01 7^ first appears. Using eqs. (|^, (|13]), one can easily check that the 
mass of the scalar field at T = Ti, = is given by 



This yields 



< - - ^ ■ (48) 



For the Higgs boson with niH ~ 60 GeV one obtains 



< >.<.^^^ - ? • (49) 



max 



Thus, even with account taken of the factor 2/3 in the expression for E, 
the dispersion of long-wave fluctuations of the scalar field is much smaller 
than the distance between the two minima. Therefore, the field on a scale 
equal to its correlation length ~ is not equally distributed between the 
two minima of the effective potential. It just fluctuates with a very small 
amplitude near the point = 0. The fraction of the volume of the universe 
filled by the field 0i due to these fluctuations (i.e. due to subcritical bubbles) 
for niH ~ 60 GeV is negligible. 

Since we already successfully applied this method for investigation of tunnel- 
ing, we expect that this estimate is also reliable. The answer remains rather 
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small even for ttih ~ 100 GeV, when the phase transition is very weakly first 
order. 



Moreover, even these long-wave fluctuations do not lead to formation 
of stable domains of space filled with the field 7^ 0, until the tempera- 
ture is below Tc and critical bubbles appear. One expects a typical sub- 
critical bubble to collapse in a time r ~ ^mix! ^^i^ is about thirteen or- 
ders of magnitude smaller than the total duration of the phase transition, 
At ~ 10~^H^^ ~ 10~^MpT~^. We do not see any mechanism which might 
increase r by such a large factor; effects such as decrease of the speed of the 
collapse of such bubbles due to finite-temperature effects (considered in the 
next section), or the inefficiency of radiating away the energy of oscillating 
subcritical bubbles are much more modest. 



Despite all these comments, we think that subcritical bubbles deserve 
further investigation. They may lead to interesting effects during phase tran- 
sitions in GUTs, since the difference between T^^ and the duration of the 
GUT phase transitions is not as great as in the electroweak case. They 
may play an important role in the description of the electroweak phase tran- 
sition as well, in models where the phase transition occurs during a time 
not much longer than T~^. This may prove to be the case for very weakly 
first order phase transitions with tuh ^ 10^ GeV, when the distance be- 
tween the two minima of V{(f),T) at T ~ Ti is smaller than the dispersion 



6 Propagation of the Bubble Wall 

6.1 General Considerations 

After bubbles are formed, one expects that they will grow until they col- 
lide. Since the expansion rate is so small at this time, provided the bubbles 
have a velocity which is not extremely small, typical bubbles will grow to a 
macroscopic size. Thus, it is important to understand how bubbles propa- 
gate long after their formation. An underlying assumption in discussions of 
the evolution of the bubble wall to date is that some time after the forma- 
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tion of the bubble, a steady state situation is achieved, in which the scalar 
field, temperature, and particle velocities are all constant in time in a frame 
which we will call the 'wall frame.' While it is plausible that such a state is 
achieved, we will not prove that this is the case; indeed, as we explain below, 
one can imagine situations in which this is not the case. To develop some 
intuition, we will consider two simple models. One picture, which suggests 
that a steady state will arise, is the following. As the wall passes through 
the medium, it dissipates energy and heats the gas. Since the wall is quite 
permeable, especially to light quarks, heat is readily transported both behind 
and in front of the wall. Once the wall is very large, the problem becomes 
one dimensional, so consider as a model a point source of heat, moving in an 
incompressible fiuid with velocity v. The temperature will obey a diffusion 
equation, with diffusion coefficient x- The solution of this equation is 



where Too is the temperature as x — > oo. Note that this is a function only of 
X — vt, so that in the wall frame one has a steady distribution, x is of order 
a mean free path, £. Thus, if v is small, the temperature varies on a scale of 
order i/v, i.e. on a scale large compared to a mean free path, i. 

Thus, it seems plausible that a steady state situation is achieved. How- 
ever, L. Susskind has given a simple example in which this is not the case, 
which we refer to as a 'snowplow.' Suppose one has, in the original, unbroken 
phase, a non-zero density of some exactly conserved quantum number; we 
will refer to these particles as 'molecules.' Suppose also, that the molecules 
have a big potential energy in the broken phase. Then it is easy to convince 
oneself that there is no steady solution; necessarily, there must be a build up 
of molecules as time evolves. Indeed, one might guess that in this case there 
is a layer of molecules in front of the wall, which grows in size linearly with 
time. Moreover, viewed in the frame of the wall, there is a steady 'wind' of 
particles. These particles collide with the ever growing layer of molecules, 
which is stationary in the frame of the wall. The molecules in this wind, 
presumably, equilibrate, so the wind comes to rest in the wall frame. This 
buildup of particles leads to an increase in pressure near the wall of order pv"^, 
where p is the molecule density and v the velocity. Balancing the stresses on 




(51) 
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the wall then gives 



v^^^j. (52) 

Note that it is important, in determining the force on the wall, that the 
properties of the gas near the wall depend on the velocity of the wall. 

In the Standard Model, there is no exactly conserved quantum number 
of this type. However, there are several approximately conserved quantum 
numbers, and one can ask whether something like this snowplow effect can 
occur. Again to develop some intuition, it is helpful to consider a simple 
model. Imagine a wall passing through a region of 'sticky dust.' The dust 
particles, when struck by the wall, stick to it, but after they hit the wall they 
have a lifetime r. Then the number of particles per unit area on the wall, 
A^, satisfies 

^ = nov--N, (53) 
at T 

where no is the density of particles and v is the wall velocity. For a steady 
state, = rioVT. In our case, the analog of the dust particles are top quarks, 
and the lifetime r is some characteristic time to change the number of tops 
and antitops (we will discuss this in more detail below). Thus, we do in 
fact expect some velocity-dependent enhancement of the top quark density 
near the wall, but we do not expect that the layer out front should grow 
indefinitely. For the rest of this paper, we will assume that the wall does 
indeed achieve a steady state. On the other hand, as we will discuss shortly, 
the enhancement of the density we have discussed above may be an important 
factor determining the wall velocity. 



Turok has argued that reflection of particles from the wall does not 
slow the wall's motion. His argument starts with the (correct) observation 
that if the velocity and temperature of the g constant across the wall, 
and if everywhere the system is described by equilibrium distributions ap- 
propriate to the local value of (f), the force on the wall is independent of the 
velocity. It is instructive to consider this part of the argument, and indeed 
a modified version of the method for estimating the force on the wall will be 
useful in our analysis. 
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In the wall frame, the net force on the gas is simply 



d^x —— — r- n(fc, x) -—— . (54) 



dt J 2A;0(27r3) ^ ' ' dx 

With the assumption that the velocity and temperature of the gas are nearly 
uniform across the wall, the density in the wall frame is simply the Lorentz- 
boosted distribution from the gas frame. On the other hand, n{k,x) trans- 
forms as a Lorentz scalar. Since 2kO(2TT'-^) ^^'^ Lorentz-invariant volume 
element, a simple change of variables gives for the net force a result indepen- 
dent of the wall velocity. Prom this Turok concludes that one needs to find 
other sources of dissipation if the wall is not to accelerate continuously. 

Our earlier discussion suggests, however, that the particle distributions 
near the wall will exhibit a more complex dependence on the wall velocity. In 
other words, even for small velocity, there will be departures from equilibrium 
proportional to the velocity. In the following sections, wc will try to take this 
velocity-dependence into account, and estimate the damping of the wall's 
motion due to scattering of particles from the wall. There will, of course, be 
other sources of damping; in some regimes, these may be larger. However, 
we will see that from this source alone, damping is sufficient to give non- 
relativistic motion of the wall for a wide range of parameters. 

The calculations which follow make use of the infrared improved effective 
potential described in Section 4 for which, as already pointed out, baryogene- 
sis has been ruled out experimentally. However, the study of the propagation 
of the bubble wall in a hot plasma receives its main motivation from the un- 
derstanding of the baryogenesis at the weak scale. We will consequently use 
the above potential as a toy model in a range of parameters so chosen that 
baryogenesis is a viable phenomenon. To this end, we will often use in what 
follows the value mn — 35 GeV for which (f)/T 1. 



6.2 Three Limiting Cases 

This discussion suggests that there are three hmiting cases which one might 
wish to consider. In the first, the wall is thin compared to mean free paths for 
all scattering processes, both elastic and inelastic. Then a typical particle. 
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as it passes through the wall, loses momentum chiefly through its interaction 
with the wall. Near the wall, one expects significant, f -dependent departures 
from equilibrium. A second extreme will occur if the mean free paths for 
both elastic and inelastic processes are short compared to the size of the 
wall. In that case, the system will be quite close to equihbrium; deviations 
from equilibrium will be proportional to the velocity and some power of a 
typical mean free path. Finally - and this is the situation which we will see is 
closest to reality - mean free paths for elastic processes may be short, but for 
inelastic processes long. In this case, in addition to the small deviations from 
equilibrium we have just mentioned, there will also be a density enhancement 
as in our snowplow discussion. We expect that the wall will be slowest if the 
'thin wall' picture is valid, fastest in the 'very thick wall' case (where mean 
free paths for inelastic processes are short) and will have some intermediate 
velocity in the third case. 

In this and the following two sections, we will consider these three cases. 
After estimating the relevant length scales, we will make a series of (admit- 
tedly preliminary) computations of the velocity and wall thickness. We will 
see that all three limits suggest a non-relativistic value of the wall velocity. 

Before considering the properties of the bubble wall, it is useful to consider 
the system at the critical temperature, Tc. At this temperature^ it should be 
possible for the phases with broken and unbroken gauge symmetry to coexist. 
Regions of different phases should be separated by a wall, which we refer to 
as a bubble wall or domain wall, at rest. It is easy to determine the form 
of this bubble wall. One is looking for a static solution of the equations of 
motion in the presence of the gas; the relevant equation is simply 

^ dV{<j>,T) 
dx^ dcj) 

We require that tend to a constant as x ^ ±cxd, so we can immediately 
solve the equation by quadrature: 



Ax = / , ^ . (56) 
V2y(0,T) 



'We note that the analysis of Ref. pOj is indeed valid atT — Tc. 
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For a Higgs mass of 35 GeV, and a top quark mass of 120 GeV, one finds 
that the wall thickness is 

5-40 . (57) 

This estimate of the wall thickness will provide a useful benchmark in what 
follows. 

Indeed, it is useful to compare this number with the mean free paths 
for various processes. In considering the properties of the bubble wall, the 
relevant mean free paths are those for particles which interact with the wall, 
i.e. principally top quarks, W^s and Z's. The processes with the shortest 
mean free paths are elastic scatterings. These exhibit the characteristic sin- 
gularities of Coulomb scattering at small angles. What actually interests us, 
however, is the momentum and energy transfer in these collisions. This is a 
problem which has been extensively studied, and we can borrow the relevant 



results. For a relativistic top quark of energy E, one has ||39 

dE Stt 



dx -y"'^ • (58) 

We expect a similar formula to hold for W and Z scattering, with as replaced 
by aw- This number is to be compared with the momentum loss due to 
interaction with the wall: 

dpx 1 dm'^ 

dx 2px dx 

Assuming that the wall thickness is of order 6 ~ 40 T~^, and noting that for 
the t quark, Am^ ~ .3T^, while for the W and Z, Am^ ~ .2T^, one sees 
that the momentum loss per unit length due to scattering is in both cases 
much larger than that due to the wall. 

We can understand this result in an alternative way. The elastic scattering 
cross section diverges at small angles, in empty space. In the plasma, we 
expect that this divergence is cut off. Examining the expression for the 



gluon propagator in Ref. we see that for the Coulomb fields, this cutoff 
is the Debye mass, 

ml = (2iV, + Nf)^^ 2.5T2 , (60) 
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while for transverse gauge bosons, it is mf^/cg/SA;^. Actually, the fact that mo 
is comparable with T means that one should go back and self-consistently 
calculate the propagators (more or less as we did in our discussion of the 
Higgs potential earher). We will not do this here, but instead simply assume 
that both the transverse and longitudinal gauge boson exchanges are cut off 
at a scale of order T. So to get a crude estimate of the mean free path, 
we simply calculate the elastic scattering cross section with a gauge boson 
propagator 

One then obtains a total cross section 

'^t^^. (62) 

To estimate the mean free path, we multiply this by the flux of quarks. 
Similarly one can estimate the contribution due to scattering from gluons, 
and the mean free paths of 1^'s and Z's. Finally one obtains an estimate for 
the mean free path of order ^ ~ 4 for quarks, and £ ~ 12 for W^s and 
Z's. We will see below that the wall velocity goes as ^/i; we expect that the 
uncertainties here will not qualitatively affect those calculations. It would 
be desirable to redo this analysis including the screening self-consistently 
and using real transport equations. However, it is reassuring that they are 
consistent with the expression above for the stopping power. 

We will also need an estimate of the mean free path for processes which 
change the number of top quarks or the number of Ws and Z's. A good 
measure of this distance is provided by the lifetimes of these particles in 
the high temperature plasma. If we treat the scalar field as approximately 
constant in space, then we expect that these decay rates have essentially the 
same form as the zero temperature rates, with the zero-temperature masses 
replaced by ^-dependent masses, and including the appropriate factor for 
time dilation. Hence the ratio of widths to masses is the same as at zero 
temperature, weighted by the Lorentz factor. These give numbers of order 
1% or smaller and thus the mean free path for decays will be of order lOO's 
of T~^'s. While we have not attempted to do a complete analysis, it appears 
from an examination of many cross sections that processes which equilibrate 
the various types of particle numbers are not likely to be effective at distances 
of order S. 
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6.3 Thin Wall 



The first limiting case we will consider, in some ways conceptually the sim- 
plest, is what we refer to as the 'thin wall.' In this limit, the thickness of the 
wall, 6, is less than a typical mean free path i for relevant particles.|3 In par- 
ticular, the momentum transferred by scattering to top quarks, W^s and Z's 
(the particles which gain mass in the broken phase) as they pass through the 
wall is small compared to the change in their mass. These particles consti- 
tute roughly 20% of the plasma. The other 80%, made up of light particles, 
are irrelevant in the sense that they interact only very weakly with the wall, 
but play an active role in establishing a steady state. In this limit, one can 
compute the force on the wall semiclassically by assuming some distribution 
on either side, and then simply following trajectories of individual particles 
across the wall. In the rest frame of the wall, energy is conserved, so it is 
easy to compute the momentum transferred by individual particles to the 
wall as they pass through or are reflected back. To estimate the force on the 
wall in this case, we will take advantage of the one dimensional geometry 
of the problem, and suppose that the wall is moving in the +x direction 
with velocity v. We will assume that the wall is surrounded by a plasma 
of roughly constant temperature and velocity. Indeed, our discussion of the 
diffusion equation (^) suggests that the temperature and other quantities 
should vary on a scale i/v. For small wall velocity, this is much larger than 
the thickness of the wall. This assumption will be justified a posteriori. 

We noted in the first section of this paper that a simple estimate, eq. (|1|), 



was obtained in Ref. |2^. This estimate is valid for strongly first order phase 
transitions, when the main difference in energy density p inside and outside 
the bubble is due to heavy particles, which acquire mass m ^ T inside the 
bubble. Such particles, coming from the phase = 0, are completely reflected 
by the bubble wall. However, in the minimal electroweak theory, the masses 
of the particles are not very large compared to the temperature and one 
cannot neglect the effect of particles crossing the wall. We have presented 
full details of the calculation in the Appendix. We quote the results here. 



^''This should not be confused with the thin wall approximation during the bubble 
formation, when the thickness of the wall is assumed to be small as compared with the 
radius of the bubble. 



38 



obtained by expanding in small v. 



v8 + = 1/(0, T) - V{ni, T) 



(63) 



Here 



S = p(0, T) - p(m, T)-—j^ E n,{E) dE . 



(64) 



In the limit of large masses {m ^ T), the last term in this equation vanishes 
and the first two accounts for the energy density contrast Ap(0, T), in agree- 
ment with (|ip. To study the properties of eqs. (p3|), ( |51D further, we make 
an expansion in powers of m/T. After a little algebra, we obtain 



S = -ET(j)' + -(j)^ 
n 4 



A 



Ar + 2fi 41n— -7 

Vo 



and 



V{0,T)-V{m,T) = -\^ 
with e given by 

rp2 _ rp2 



1-el 



In the case of small e, when T is close to Tc and 6 ~ </)c(l + e) 



At At ^ Vo 




(65) 

(66) 
(67) 

eq. (|63|) gives 
(68) 



The relevant value of e corresponds to the time when the biggest bubbles 
propagate in the plasma, which is also the time when bubbles fill up the 
universe. Fig. 6 shows that e varies between 0.1 and 0.3 for a large range of 
Higgs masses and is rather insensitive to the top mass; we use the value 0.26 
corresponding to a Higgs mass of 35 GeV to illustrate our results. For a zero 
top mass,Q the term in parenthesis is approximately 1; in such a case, 



TT 

V 

6 1 



0.1 . 



(69) 



^^This is clearly an unphysical value. However, this information will be useful in the 
discussion of thick walls in the next subsection, since for thick walls the top quark contri- 
bution is relatively suppressed. 
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As the top mass increases to a value > m^, its contribution becomes dom- 
inant and the velocity decreases significantly; as an example, for mj = 120 
GeV, V ~ 0.06. We have analyzed eq. ( |63D numerically, Fig. 7 shows the 
dependence of the velocity on e while Fig. 8 illustrates the dependence of 
the velocity on the top mass. 

We now elaborate further on the situation describing the steady state. 
Since particles change their momentum in crossing or bouncing off the wall, 
the local thermal distribution is spoiled; its restoration requires the release of 
energy-momentum at a certain rate, implemented by the fiow of light quarks. 
Clearly, to satisfy the conservation of energy and momentum, the plasma has 
to adjust its velocity and temperature distributions. For this purpose, the 
wall and the non-thermalized particles can be viewed as a source of energy- 
momentum of a size of a few mean free paths, moving in a relativistic plasma 
with a non-relativistic velocity. This situation is similar to the diffusion 
problem described in Section 6.1. Qualitatively, one expects a steady state 
situation in the rest frame of the source, a rather uniform distribution of 
temperature and velocity behind the source extending across the source on 
a typical distance ~ ^. From our estimate of the velocity above. Id is 
typically 20 times larger than i which justifies the assumption that the wall is 
interacting with a plasma of uniform velocity and temperature. Furthermore, 
as a good approximation, we can assume this velocity and temperature to 
be the ones far away from the wall. We can then write equations for the 
conservation of energy and momentum across the wall: 

i\v'^p{mX) - V{mX)) - 7'(^V(0,T) - \/(0,T)) = , 

i^v'[p{m, T) - y (m, T')) - 7'^(p(0, T) - \/(0, T)) = , (70) 

where v' = v + 5v and T' = T + 6T are the quantities defined behind the 
wall. 

It is easy to compute the changes in velocity 6v and temperature 6T 
of the plasma as it crosses the wall. These equations and the steady state 
assumption, together with equation 

{V{0,T') -V{m,T')) -v'S' c:^0 , (71) 

yield a unique answer for these quantities. Typically, ^ ~ +1%, so the ve- 
locity, V, is still given accurately by Fig. 7. For the temperature variation. 
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one finds ^ is negative and of tlie order of or less than — 0.0001%.|3 Fur- 
tliermore, 5v is positive; the process is a deflagration. It is also easy to check 
that a small amount of entropy is produced in the plasma, in agreement with 
the second law of thermodynamics. 

One can attempt to obtain the shape of the wall by writing eqns. ([70|), 
( |7TD in a differential form. At e = 0, the resulting equation is just that which 
we studied earlier for the domain wall separating the coexisting phases. At 
finite e, however, there are velocity-dependent corrections to the equation, 
and it is more difficult to analyze. Because the temperature and velocity 
vary in space, the problem one has to solve is analogous to the motion of 
a particle in a time- dependent potential. Moreover, the two minima of this 
potential are not degenerate, so one needs to understand how the oscillations 
of the scalar field in the true minimum are damped. However, for small e 
(and f ), we do not expect the shape of the wall to be too much different than 
for e = 0; in particular, the size of the wall (see eq. (|56D ) is still given to a 
good approximation by 



E 



As this is a several times bigger than the mean free paths of the relevant 
particles {('.w,z ~ 12 and Itop ~ 4T~^), the thin wall approximation 
just described is questionable within the minimal version of the electroweak 
theory. However, this approximation may be valid in some extension of the 



standard model. For instance, the authors of Ref. |]22| have argued that 
baryogenesis at the weak scale is viable in a model with a singlet scalar and 
one Higgs doublet. The presence of a light singlet makes the transition more 
strongly first order, and, consequently, allows for a Higgs mass above the 
experimental bound. In such models, the wall tends to be thinner. 

So far, we have assumed for our thin wall analysis that the density to 
the left and right of the wall (of particles moving to the right and left, re- 
spectively) are precisely the equilibrium densities. However, our discussion 
earlier of the snowplow problem suggests that there may be other effects we 
should consider as well. Suppose that there were only elastic processes, i.e. 



^^The gradient of temperature, goes like ^ ~ —v5v. 



41 



no processes which changed the separate numbers of different particles and 
antiparticles. Suppose, also, that (for example) all top quarks were reflected 
from the wall. Then we would be in precisely the snowplow situation, except 
now in a relativistic version. 

More realistically, consider the possibility that, say for top quarks, the 
probability of reflection from the wall of a given quark is /. Suppose, also, 
that the mean free path for processes which change the number of tops or 
anti-tops is r. Then there will be some buildup of top quarks near the wall. 
The time required for the wall to catch up with a given top quark goes as 

where / is a typical mean free path. As a result, in the limit of very small 
V, the wall doesn't catch up with a typical particle before it undergoes an 
"identity change." Thus, we have something like the "sticky dust" picture 
described earlier. Per unit area of the wall, there is an increase in the number 
of top quarks of order rfvn, where n is the equilibrium top quark density. 
These quarks are spread over a distance of order \/£r, so their density is of 
order fvn^. This gives an extra contribution to the force on the wall of 

order f Apn/y^, where Ap is the free energy difference on the two sides of 

the wall. It is easy to compute /; for bosons, / = for fermions, the 

result is 2/3 as large. Assuming that the square root in this expression is 
a number of order 5-10, this is comparable to the force we have computed 
above, and will tend to decrease the velocity of the wall. This discussion, of 
course, is extremely crude, but it suggests that there are various effects at 
least as large as those we have considered above, all of which slow the wall. 
Thus, if the thin wall approximation is a good guide, the wall is likely to be 
quite non-relativistic. We will see that there are similar effects in the the 
case of a thicker wall. 



6.4 Thick Wall 

We now consider the case that the wall is extremely thick. As is clear from 
our earlier discussion of mean free paths, this, like the thin wall case, is not 
completely realistic. These calculations, however, should bracket the true 
situation. At the end of this section we will try to estimate the effects of 
the density enhancement which occurs because particle numbers are approx- 
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imately conserved. So we first consider the case where the mean free paths 
for both elastic and inelastic processes are short. To get an idea of how finite 
elastic scattering lengths affect the velocity of the bubble wall, we assume 
that particles propagate freely over distances of order a mean free path, I. 
We view the bubble wall as a succession of slices with thickness of order ^, 
and for each of these we repeat the thin wall analysis. We refer the reader to 
the Appendix for the details of the derivation and the precise formulae, and 
just summarize the results here. We write the result in the following form 

e = 56£:f ^" + SfEf''' , (73) 

where <S are suppression factors dependent on I. They are not very sensitive 
to rrit and vriH-, and behave as {i/Sy^'^ (see Fig. 9). Using the values quoted 
above {f.w,z ~ l2T-\ itop - 4T-i and S - 40 T-^), we find Sb ~ 0.45 and 
Sf ~ 0.15. Using the equations for S from the previous section and assuming 
niH ~ 35 GeV and rrit ~ 120 GeV, we find a velocity of about 

~ 0.2 . (74) 

Using the values above for Sb and Sf, Fig. 10 illustrates the velocity for a 
range of top masses. 

In reality, however, we expect, because the numbers of tops and anti-tops 
(and similarly W^'s and W^'s) are not quickly equilibrated, their densities 
will be enhanced near the wall; since this enhancement will depend on the 
velocity, there will be a velocity-dependent drag on the wall, similar to that 
we discussed above in the thin wall case. 

Again, we will content ourselves with an extremely crude estimate. As 
in the thin wall case, we imagine that as particles pass through the wall, 
some fraction is refiected. These refiected particles make a random walk 
until they decay, or undergo scattering processes which change their identity, 
with a lifetime r. Then the number of particles per unit area on the wall is 
N — fnvr, where as before, / is the fraction of particles which are reflected. 
This extra density of particles will be spread over the thickness of the wall, 
or over ^/Jt, whichever is larger. For our crude estimate, we will assume this 
density is spread uniformly over the wall. In other words, we will assume that 
the density is the equilibrium density (for a given value of the scalar field) , 
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except enhanced by a factor of the form 1 + In this case, one obtains 
an increase force on the wall of order AF = ^^y^ x Ap, where, as before, Ap 
denotes the internal energy difference on the two sides of the wall. As an 
estimate of the quantity /, we take the ratio of the equilibrium densities on 
the two sides of the wall. For top quarks, this gives a number of order 5%. 
The ratio j is likely to be of order 5 — 10. Thus, this effect may be even more 
important than the effects we have considered above, i.e. we have probably 
overestimated the wall velocity. 

We can also ask about the wall thickness in this limit. Our discussion 
here suggests that the thickness will be modified from its value at To by a 
factor of the form 1 + j g i^y ^ small amount. However, it is somewhat 
harder to develop a complete theory of the wall shape in this hmit. For 
example, in the analog mechanics problem, if one allows for spatial variation 
of the temperature, one has to consider a system with time-dependent forces. 
Moreover, one has to consider how the 'motion' of damps out in the region 
of broken symmetry. This may require consideration of types of damping 
other than those we have considered up to now. We will not explore this 
issue further here, and simply assume that the shape of the wall is only 
slightly modified from its form at Tg. 



7 Conclusions 



The study of the electroweak phase transition began two decades ago, and 
stimulated work in many areas of what has come to be called astroparticle 
physics. However, until now it was not very important to know any details 
of the theory of this phase transition. For most applications it was quite 
sufficient to know that in the early universe at a temperature higher than 
about 10^ GeV the symmetry between weak and electromagnetic interactions 
was restored. Recently it has become clear that if we wish to investigate the 
possibility of electroweak baryogenesis, we must have a complete and detailed 
picture of the phase transition, from the accurate computation of the critical 
temperature to the investigation of the motion of the bubble walls. In this 
paper, we have taken some steps towards a systematic investigation of all 
relevant features of the first order phase transitions in electroweak theory. 
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We have seen how to organize the perturbation expansion in hght of the 
infrared problems which exist at high temperature. We have shown that no 
dangerous hnear terms arise in the effective potential; on the other hand, we 
have seen that the coefficient of the cubic term, which is responsible for the 
first order nature of the transition, is reduced by a factor 2/3 from its lowest 
order value. This means, in particular, that theories with a single Higgs 
doublet cannot be responsible for the observed baryon asymmetry, given the 
present experimental limits on the Higgs mass. 

We have also understood some aspects of the strongly first order phase 
transition relevant to baryogenesis. For such theories, we have seen that the 
phase transition typically proceeds through the formation of critical bubbles 
with thick walls. We have developed a method of analytic investigation of the 
probability of bubble formation, which is valid for a large class of theories. 
With the help of the stochastic approach to tunneling, we have found that 
subcritical bubbles are only likely to be important (in the minimal standard 
model) for Higgs mass larger than about 100 GeV, when the phase transition 
is second order or very weakly first order. 

We have considered the problem of propagation of the bubble wall. To 
this end, we have considered the minimal Higgs model with a light Higgs 
{niH < 35 GeV). While this possibility is ruled out, we believe this theory 
is a good toy model, whose features mimic those of more realistic theories 
in which baryogenesis is possible. Investigation of the bubble wall motion 
turns out to be surprisingly difficult; indeed, we have identified important 
mechanisms for slowing the wall which have been omitted from previous 
treatments. Our own treatment is crude, and may omit additional damping 
effects, but it suggests that the bubble wall is typically non-relativistic or 
only mildly relativistic in these theories, and tends to be rather thick. It 
would be interesting to construct theories where this is not the case, which 
might realize the scenario of ref. IT^ in which baryons are produced quite 



efficiently. 

While in this paper we have seen that theories with a single Higgs doublet 
cannot give rise to the observed asymmetry, there is still much work to be 
done to determine whether or not baryogenesis can occur in extensions of the 
minimal theory. Needless to say, such extensions are interesting in their own 
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right, and also because they can provide larger CP violation than exists in the 
minimal model. Apart from the issues of the phase transition discussed here, 
further work on B-violation rates and the detailed mechanisms of baryon 
production is still necessary. Hopefully the observations contained in this 
paper will represent a positive step on the road to a complete understanding. 
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8 Appendix: Calculations of the Wall Veloc- 
ity 



In this appendix we give the details of our calculation of the force exerted on 
the advancing bubble wall by the plasma. As discussed in the text, we will 
assume that there exists a steady state; that is, there is a well-defined rest 
frame of the wall, at sufficiently late times after the appearance of the bub- 
ble. In this frame, we have a (time-independent) velocity and temperature 
distribution for each component of the plasma that reduce to those of the 
surrounding Universe at spatial infinity. We focus on those species that are 
relatively heavy in the broken phase; these will obviously give the greatest 
contribution.^ 

In principle to solve this problem in all generality, we would write Boltz- 
mann equations for each component of the plasma, as well as an equation of 
motion for the scalar fields. However, in the simplest cases, we can learn a 
great deal by looking at the equations for local energy-momentum conserva- 
tion. These may be written in the form 

d.T^" + d.T^Z = , (A.l) 

where we have (arbitrarily) separated off the classical (zero-temperature) 
stress tensor of the scalar field. We can also write 

- Fi^) = 9.T(f , (A.2) 

where (i) labels components of the plasma and Fj*^ is the force density on the 
wall of the i^^ species. Combining ( |A.1| ) and ( |A.2D , we find the full equation 
for the scalar field: 

W = E f!^^ ■ (A.3) 

i 

We now turn to the evaluation of F^.. 

Consider a volume element of width s in the x-direction. For our analysis, 
we will assume that particles of the plasma typically traverse this distance 



^•^The light degrees of freedom can not be ignored however, as they are instrumental in 
the approach to and establishment of local kinetic equilibrium. 
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without interacting; we can calculate the force on this slice of the bubble 
wall by following individual particles of the distribution. We will present the 
analysis for a general s below, but note now that the size of s distinguishes the 
three cases considered in the text. Namely, the thin wall scenario corresponds 
to s --^ S. where S is the wall thickness, the very thick wall corresponds to s 
infinitesimal, and the intermediate scenario to s ~ with i some relevant 
mean free path. We will suppose that particles entering the volume element 
on either side are described by equilibrium distributions appropriate to the 
masses on either side {nio or mi, with mi > mo). Note that particles leaving 
the volume element will not, in general, be thermally distributed; we assume 
that equilibrium is restored in a distance of order a mean free path. In other 
words, for particles entering the volume, we are assuming the distribution 

n{E) = no{j,{E - vpx),T, fj) . (A.4) 

with E — yjp^ + mil- 

We compute the force per unit area as follows. In time dt, there are 
n{p, T)^;^\vx\dt particles per unit area providing a force Ap^/dt. Thus, the 
pressure is 

dP^I ^^n{p,T){Apx)\vx\ . (A.5) 

We write the integration measure as dE E dpx/4:7i^. There are various re- 
gions of integration, corresponding to whether the particles have sufficient 
momentum {pi > {ml — ml)) to go through the slice, or are reflected. 

For particles bouncing off the wall, incident from the right, the momentum 
transfer is just —2px. We arrive at 

APi^^J dE dpx plno{jv{E - vpx),T) , (A.6) 



where the region of integration is [{mo, mi) x (— w (E'^ — ml), 0)] + [(mi, oo) x 



(— cr, 0)] with a = yml — m^. 

For particles incident from the right with \px\ > o, we flnd Ap^ — 
~\JPx ~ ~ Px a-iid the pressure is 

APii = j ^dpx Px (^Px + \/pI- 0-^) no(yv{E - vpx), t) , (A.7) 
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where the region of integration is [(mi, oo) x ( — y (E"^ — ml), —a)]. 

For particles incident from the left, the momentum transfer is Ap^ ■ 
pI + (j"^ — Px and so 

APiii = ^JdE dp, px (^-p, + ^pI + cr^^ no{-f,{E - vp,),T) , (Ai 



where here the integration region is [(mi, oo) x (0, ^{E'^ — mf))]. 

We need to now integrate over the thickness of the wall. If there are 
N 6/s slices, we write 

TV 

Fx = J2 ^Pn , (A.9) 

ra=0 

where here AP„ is the sum of eqs. ( |A.6| )-( |X^ ) with m^, the mass at the rfi^ 
step, given by 

rrio = dx——— ~ m — , (A. 10) 



oo 



dx m 

In the last step we have made a linear approximation to the wall profile. 

We may proceed in a variety of ways. If we assume that the speed of the 
wall is small, we can expand the distribution ( |A.4| ) to lowest order in v. On 
the other hand, we could expand in powers of m/T, at the risk of making a 
small error for the top quark. In this Appendix, we will choose the former 
expansion. We find 



n, 



o{i.{E - vpx)) = no{E) - vpx ^{E) + 0{v') . (A.ll) 



To lowest order we have: 

Ai^°l = \r^,n^{E,T){E'-mlf'^-{m^^m,) 

6 Jnio ^TX 

= Fimo,T) - F{mi,T) . (A.12) 
Summing over the wall and over species, we find the simple result 

F]o] = F(0) - F(m') . (A. 13) 
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The pressure difference, to lowest order, is the difference of free energy den- 
sities. We have neglected contributions proportional to gradients in temper- 
ature and chemical potential. 

The next-to-leading order terms account for the forces of the plasma on 
the wall. To order v^, we find 



-2v 



^1 dE 



—V 



mo 47r2 
~ dE 



noiE)E{E^ -m[ 



(A. 14) 



.1 47r2 



n,(E) E (^^ (E^ ~ ml) - ^ {E^ - ml)^ . 



The next step is to sum this over the wall. We will do this for the three cases 
separately. For the thin wall, we simply take = 1, rrio = and mi = m, 
and obtain 



= V ip{0)-p{m) -m^ 



°° dE 



47r2 



Eno{E) = v£ 



(A.15) 



It is useful to expand this expression in powers of m/T. For W and Z bosons, 
we find 



TT 



647r2t;4 



2mt/ In 



+ In 



Likewise, for top quarks, we findj^ 



aBvlT^ 2^ 

(A.16) 



£ 



thin 



W 4 



^/ levrV""* 
The total £ is then 



In 



7 



apv^T^ 2 



+ 0((^)^) 



(A.17) 



£ ~ -ET(j)^ + -0^ 

TT 4 



A - At + 25 4 In 



(A.18) 



If the wall is much thicker than any relevant mean free path, we take s 
proportional to dx. We can easily see that (A. 14) has no linear term in dx, 

^*The expansion is less trustworthy here; in particular, we have neglected terms of order 
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and we recover Turok's result that there is no w-dependence in the force 
on the wall in the limit £ — 0.0 Presumably there are other effects that are 
ignored in this computation that slow down the wall. 

The case of greatest interest is finite Ijb. We have performed the analysis 
numerically and find E is suppressed essentially by a factor of (f/^)^/^. (This 
non-analytic behavior is expected of eq. (A. 14), since its Taylor expansion 
has singularities at order cx^.) We can now write E in the following form: 

E = SfEf'' + f . (A. 19) 

where 

- e 



Sf and Sb are rather insensitive to rrit and ttlh- We have plotted these factors 
in Fig. 9 vs. {i/Sy^"^. The behavior of this plot is adequately approximated 
by a linear function: 

S^ - [-^j (A.21) 

for both fermions and bosons. We find then, that the wall velocity, for 
niH = 35 GeV and rrit = 120 GeV, is approximately v ~ 0.2 (see Fig. 10). 



^^We have explicitly checked this through order w^. 
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